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Q> ' transition b — > S77, including leading order QCD effects. We apply the dis- 

Q_i' cussion to the decay Bg 77. Even though one-particle irreducible con- 



We study the effect of supersymmetric contributions to the effective quark 



tributions could play a role, numerical cancelations make the amplitude for 
the two-photon emission strongly correlated to the 6 — > 57 amplitude which 



^ ■ is sharply constrained by experiment. A quite general statement follows: as 



long as non standard physics effects appear only in the matching of the Wil- 
son coefficients of the standard effective operator basis, the deviations from 
the standard model expectations of the decay rates induced by 6 — > 577 are 
bound to follow closely the corresponding deviations on 6 — > 57. Effects of 
new physics are therefore bound to be small. 
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I. INTRODUCTION 



The rare flavour changing transition b s'~)'~) has recently attracted new interest in 
view of the planned experiments at the upcoming KEK and SLAC B-factories and existing 
hadronic accelerators, which may test branching fractions as low as 10~^ times the B meson 
decay width. 

Rare B-physics potentially provides us with valuable redundacy tests of the flavour struc- 
ture of the standard model (SM) and complementary information on the related CP viola- 
tion. It is also most sensitive to the "heavy" sector of the SM particle spectrum and it is 
the preferred low-energy laboratory for virtual signals of exotic physics. 

In recent years the 6 — 57 decay has provided us with the first experimental evidence 
of "penguin" physics and the experimental bounds on the B —>■ Xg'j decays have shown to 
provide sharp constraints on the physics beyond the SM. 

The fact that the 6 — S7 transition has been already experimentally observed is related 
to the peculiar enhancement of the electroweak amplitude that arises at the two-loop level 
due to a large logarithmic QCD mixing with the effective b sec operator [0,0 . The study 
of the QCD leading logarithmic (LO) resummation for b ^ s'y (and 6 — > s gluon) has 
been very recently extended to the next-to-leading order thus reducing the theoretical 
estimated error for the inclusive rate below the 10% threshold. 

In order to produce similar studies for the b 577 transition it is extremely helpful to 
observe that, by the use of an extension of Low's low energy theorem |^,^ or, alternatively, 
by applying the equation of motions, the b 577 quark operator can be expanded at 0{G f) 
on the standard operator basis needed for b — > 37. 

Three groups have recently presented a QCD LO calculation of the two photon transi- 
tion B^10| thus improving on the previous electroweak calculations [^|1T],|12] . 



More recently a study of the b 577 decay in two-Higgs doublet extensions of the SM 
has appeared . 



The purpose of the present paper is to study the influence on the b —>■ S77 transition 
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of softly broken supersymmetry. We find actually the results of our analysis having a more 
general impact on the possible effects of new physics for the radiative two-photon decay. At 
the LO the short-distance part of the b S77 amplitude turns out to be controlled by the 
one-photon radiative component. The higher dimension one-particle irreducible contribu- 
tions present in the h 577 amplitude, potentially large because of the I /ml dependence, 
turn out to remain subleading because of accidental cancelations. This result remains true 
when studying two-Higgs doublet extensions of the SM where the additional charged Higgs 
contribution to the 6 ^ 57 component of the two-photon amplitude adds coherently to 
the SM one. On the other hand, in supersymmetric theories there are potentially large 
destructive effects related to the exchange of superpartners, which may reduce the size of 
the one-particle reducible part of the b 577 amplitude thus affecting the relative weight 
between the latter and the one-particle irreducible components. 

Nevertheless, the present experimental constraints on 6 ^ 57 induced decays are tight, 
and by studying the effects of low energy supersymmetry as a paradigmatic case, the present 
analysis shows generally that, as long as extensions of the SM affect only the short distance 
Wilson coefficients of the standard effective Hamiltonian, the b 577 induced decays are 
severely constrained by the present bounds on the inclusive B X^j decay 

BR{B X,7) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10"^ (1) 

Since this result is consistent within 30% with the next-to-leading order (NLO) SM calcu- 
lation |] 

BR{B X,7) = (3.28 ± 0.33) x 10-^ (2) 

one may not expect much larger deviations of the two-photon decay rates from the cor- 
responding SM estimates. Were one to invoke supersymmetry to reduce by 30% the SM 
prediction of r(i? —>■ Xg'-f) then the SUSY rates related to 6 ^ S77 will unlikely exceed 
the corresponding SM expectations. Potential implications of a (albeit challenging) NLO 
analysis of the two-photon amplitude in extensions of the SM are commented upon in the 
conclusions. 



Although we shall apply our results to the Bg 77 decay we are not interested in 
predicting absolute decay rates and thus we will not discuss the uncertainties related to 
hadronization, which are not affected by short-distance new physics. Our main purpose is 
to study the deviations of the quark amplitudes from the SM predictions. 

On the basis of these considerations and at present time, it is not worthy the effort to 
perform a detailed analysis of a specific supersymmetric model. We will investigate the 
main features of the SUSY amplitudes by means of relatively simple limiting cases in the 
supersymmetric and soft breaking parameter space. 

II. EFFECTIVE QUARK LAGRANGIAN: OPERATOR BASIS AND 

COEFFICIENTS 

At the LO in QCD the effective Hamiltonian for 6 ^ 57 closes on a basis of eight effective 
operators 

Gp ^ 



H, 



i=l 



(3) 



where 
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Sia'^"(m,(l - 75) + mb{l + -ir>))biF^ 



- 75) + m,{l + i,))^b,G% 



(4) 



In eq. (|) i,j = 1,2,3 are color indices, a = 1,...,8 labels the SU{3) generators, and 
± A = 1 ± 75. The sum runs over the active quark flavours u, d, s, c, b. 



Having factored out the relevant Kobayashi-Maskawa (KM) mixings, the LO matching 
of the Wilson coefficients Ci(/x) at the scale mw is given in the SM by 



C^^\mw) = 0, i = 1,3,4,5,6 (5) 



Cr{mw) = l (6) 
('"»■' = - 24(.-l)3 P' 

where x = m^/m^. 

Large (a^ log /i)" corrections to the weak amplitudes are resummed via the renormaliza- 
tion group (RG) equations by evolving the Wilson coefficients at the typical scale of the 
process (/i ~ mi,). While the anomalous dimension matrix of the first six operators is reg- 
ularization scheme independent, the operator mixings between Oi, ■ ■ ■ , Oq, and the dipole 
penguins Oy^s are generally regularization-scheme dependent since they arise in LO at the 
two-loop level. However, the finite one-loop matrix elements of 6 — > 57 (and b sg) gener- 
ated by the insertion of the O^^ gluonic penguins carry a compensating scheme-dependence 
such that the total physical amplitudes are independent on the regularization scheme. In 
view of this, one defines the so called "effective" Wilson coefficients p^JT6|] , C7-^-^(/i) and 
Cl^^ifi), for which the LO RG running is scheme independent. As an example, for Crlfi) 
one finds the LO expression 

16 8 14 16 ^ 

Cr^M = V^Crimw) + - (r/^ - r?^! ) Csiniw) + C2{mw) E (9) 

^ 1=1 

where t] = as{mw)/c(s{f^) and the numbers hi and are given in the appendix of ref. JIB 



The sum of all hi is zero and Ci^^{mw) = Ci{mw)- For notational convenience, the super- 
script "eff" on Cj^sifJ') "will be henceforth omitted. 
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III. THE Bs 77 DECAY 



The quark b 577 transition induces at the hadronic level two interesting rare decay 
modes of the B mesons: B^ ^i Xg 77 and B^ 77, where Xg represents strange mesonic 
states. Both these decay modes and their features at the hadronic level within the SM have 
been widely studied in the literature [P!7| . 

We apply the analysis of the QCD corrected supersymmetric b —>■ S77 amplitude to 
the discussion of the two body Bg —>■ 77 decay, whose total rate can be cast in a simple 
and compact form; this allows us to carry out the present analysis without unnecessary 
complications. Our conclusions are based on short-distance properties of the quark transition 
that hold as well for the Bu,d Xs 77 decay, even though the detailed form of the amplitude 
differs from that of the Bg 77 decay. 

The total Bg — > 77 amplitude can be separated into a CP-even and a CP-odd part 

A{Bg ^ 77) = M+F^^F'^'' + iM-F^^F'^r (10) 
According to the notation of ref. one finds 

= {b{,) K{ml) + ^) , (11) 

and 

M- = ^^^^^^^fB.msy,:V,, (j:mn^ A,(/i) J(mJ) + B{^,) L{ml) + , (12) 

where A = niB— and 



An = (C3 - C,)N, + (C4 - Q 
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A, = ^iiC,-C,)N, + {C,-Ce)) 
A, = (Ci + C3 - C,)N, + (C2 + C,-Q 
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A = A, = ^ ((C3 + C4 - C5)iVe + (C3 + C4 - Ce)) 
B = -]{CeN, + C5) , (13) 



are combinations of Wilson coefficients evaluated at the scale ^ (~ rrib). 
The functions J{m?')^ K{m'^) and L[m?') are defined by 

J(m2) = /ii(m2), 
ir(m=^)=4/n(m2)-/oo(m2) , 

L{m^) = looim'), (14) 



with 



Ipq{m^) = dx dy— — , (15) 

JO Jo m'^ — 2xyki ■ k2 — te 



,2 



and 2 ki ■ k2 = ms^- 

The decay width for Bg —>■ '~f'~f is then given by 

T{B, 77) = ^(|M+|' + \M-f) , (16) 

and, from the measured Bg lifetime, the corresponding branching ratio is finally obtained 

BR{Bg ^ 77) = r{Bg ^ 77)/r(5,) . (17) 

Coming to the inclusive Bg — > X^'j decay it is convenient to use the approximate equality 

BRiB ^ X,7) ~ ^}!''"'^_\ bR{B ^ X,ez>e) , (18) 
r(o — > cez/e) 

where 

rfh ^ cr.A IT/*T/:.|2 

ICri^^W , (19) 





si) 


_ \V*Vtb\^ 6aem 


T{b - 


ceue) 


iVcbP 7ig{z) 



which minimizes the uncertainties related to the bottom quark mass and KM mixings. In 
eq. (|19|) the function g{z) = 1 — + 8z^ — — 2Az^ log z is the phase space factor in the 
semileptonic decay and z = rric/mb- 

In Figs. 1 and 2 we show the LO results for the SM Bg 77 and B —>■ Xsj branching 
ratios, as a function of the renormalization scale /i and of rrit. Our numerical results are 
obtained using the values given in Table | for the other input parameters. 



Fig. 1. BR{Bs ll)sM X 10^ as a function of /z and rrit (GeV) for central values of the 
other input parameters (Table I). 
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BR( b -> sY ) 2.5 

2.25 




Fig. 2. BR{B Xs'y)sM x 10"^ as a function of /i and (GeV) for central values of 
the other input parameters (Table I). 

From a direct comparison of Figs. 1 and 2 it appears clear that the Bg 77 decay rate is 
dominated by the Cj component. This is related to the fact that the one-particle irreducible 
contributions arising from the operators Qi^2, which could be potentially large due to the 
1/ml dependence of eq. (|T5|), appear in eq. (|13D (Ac) via the combination Ci + C2 which 
is numerically suppressed H. Notice that the scale dependence represents the largest source 



of uncertainty of the LO calculation ||16| , p!8| . We have here shown the range 2.5 GeV < < 
10 GeV. 

As a reference for the following analysis, the SM central values for the LO QCD corrected 
decay rates at the scale /i = mj, are given by 

BR{B ^ Xs-f)sM = 2.5 X 10-" (20) 
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and 

BR{Bs ^ i-i)sM = 4.4 X 10-^ . (21) 
The prediction in eq. (pT]) compares to the present experimental bound |19 



BR{Bs 77) < 1.48 X 10"^ (22) 
which is about three orders of magnitude away from the needed sensitivity. 

IV. THE Bs 77 DECAY IN SOFTLY BROKEN SUPERSYMMETRY 

In a wide class of realistic SUSY models the global supersymmetry breaking is a conse- 
quence of the spontaneous breaking of an underlying N=l supergravity theory (for reviews 



see ref. ^^). The locally supersymmetric lagrangian is supposed to undergo a spontaneous 
breaking in the so called hidden sector, and the effects of this breaking are communicated to 
the observable sector through gravitational effects. A renormalizable theory is obtained in 
the limit in which the Planck mass goes to infinity. By doing so one is left with an effective 
globally supersymmetric lagrangian and explicit soft breaking terms. In our present study 
we shall consider the following gauge invariant soft breaking Lagrangian: 

= -M^ -{M + S + h.c.) (23) 

where is a common mass term for all the scalar components Zi in the theory 

= ^im^z*Zi , (24) 

M is a mass term for the gauginos Aq, a = 1, 2, 3 considered as Weyl fields 

M ^ -^KK , (25) 

and S is the scalar analogue of the superpotential 

S = m \-AuhuH2QU^ + AohDHiQlD^ + AEhEHiLE" + B^iHiH^ , (26) 



where with standard notation hu,D,E are the 3x3 Yukawa matrices for the quarks and 
charged leptons. The soft breaking parameters and B are dimensionless numbers of 
order unity. 

The b —>■ sj and b —>■ S77 transitions can proceed in the SUSY model via five different 
intermediate particles exchanges: 

1. Charged gauge bosons iW~) + up-quarks 

2. Charged Higgs bosons {H~) + up-quarks 

3. Charginos + up-squarks 

4. Gluinos [g) + down-squarks 

5. Neutralinos (x^) + down-squarks 

The total amplitude is the sum of all these contributions. The complete analytic expres- 



sions for the various components are found in ref. pl |. 

An effective b — s flavour changing transition induced by W~ exchange is the only way 
through which the decays proceed in the SM. A two-Higgs doublet extension of the SM 
would include the first two contributions, while the last three are genuinely supersymmetric 
in nature. 

Gluinos and neutralinos can mediate flavour changing interactions only via renormal- 
ization effects which are crucially dependent on the detailed structure of the model. Their 
consideration is beyond the scope of the present work and our results do not presently jus- 
tify a more detailed analysis. We shall discuss the features of the inclusion of the first three 
contributions in the matching of the Wilson coefficients. 

The supersymmetric Wilson coefficients are then given by 

C^S'^^imw) = Cl^imw) + C^,si^w) + CUmw) , (27) 
while at the LO the matching conditions in eqs. (^-(§P remain unaffected. 
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From the results of ref. |^T| and comparing with eq. we obtain the following contri- 
butions to the Cf^si^w) coefficients: 



H 



tan^/3 



(1) 

7,8 



m 



(2) 



H , 



(28) 



induced by charged Higgs exchange and 



EE 



m 



w 



y'^UL ~ J^URJK^UL ~ J^UR )J7,i 
^2 



T*jks\ r{l) 



TTjkb/^*jks TT*jks\ '^Xj f (3) ( 



(29) 



induced by chargino exchange. We have found convenient for the present discussion to 



introduce the functions f^^^ according to the notation of ref. 



(1 



f: 



f? 



f- 



(3 



(1 



(2 



fi 



(3 



X 



X 



X 



X 



[x 



(7 - 5x - 8x2) ^(^3^ _ 2) 
36(x- 1)3 ^ 6(x-l)4 
(3 - 5x) (3x - 2) 



logx 



3rx- 1)3 



6(x- 1)2 

(2 + 5x - x2 

12(x - 1)3 

(3-x) _ 1 
2(x - 1)2 ^ (x - 1) 



logx 



X 



2fx - 11 



■ logx 



logx 



(l-x)/f (x)-^/8^^^(x)-^. 



X 



•(2)/ 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



These functions have simple and obvious relations with the functions F^ix) defined originally 



in ref. 211, to which we refer the reader for all details. 



In eq. (|29|) j = 1,2 is the label of the chargino mass eigenstates and k = 1,...,6 is 
the analogous label for the up-squarks; the matricial couplings Gul arise from charged 
gaugino-squark- quark vertices, whereas Hm and Hur are related to the charged higgsino- 
squark-quark vertices. These couplings contain among else the unitary rotations U and V 
which diagonalize the chargino mass matrix 

M2 mvi/A/2 sin/5' 



U* 



mw 



72 cos/? 



V 



-1 



-/i 



m 



(36) 



X2 



11 



where M2 is the weak gaugino mass and /i the Higgs mixing parameter. The sign of the n 
entry is defined accordingly to the Feynman rules used in obtaining the above results (see 
the comment following eq. (15) in ref. p3|). 

Due to the relevance of the chargino amplitude for the present discussion, it is worth 
trying to have a better understanding of the nature of the features exhibited by this ampli- 
tude. 

An explicit tan/3 dependence is found in Hul and Hur where quark Yukawa couplings 
are present; more precisely, Hjjl is proportional to the down-quark Yukawa coupling, which 
grows with tan/3 as 1/ cos/3, whereas Hur contains the up-quark Yukawa coupling, that 
approaches in the large tan (3 limit a constant value (oc 1/ sin (3) . It is in fact the contribution 
in the second line of eq. (^) that determines the behaviour of the amplitude in the large 
tan f3 regime. Detailed studies of this feature of the chargino amplitude are available in the 



literature 24 



In order to allow for an analytic and more transparent discussion of the chargino compo- 
nent we resort to simplified assumptions on the squark mass spectrum which reproduce with 
good approximation the global features of the model. In this we follow closely the analysis 
of ref. p2[. We assume that all squarks, other than the two scalar partners of the top quark, 
are degenerate with the soft breaking mass m. The remnant 2x2 top squark mass matrix 
is diagonalized by an orthogonal matrix T such that: 

+ {Atrh + fi/ tan (3)mt 

{Atm + fi / tan f3)mt m? + rn^ 

where At is the supersymmetry-breaking trilinear coupling. The sign of the /i term is con- 
sistent with eq. (p^). 

With these assumptions eq. (|29|) can be written as 



T 



T- 



1 _ 


(0 









ml 



(37) 
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-E 

k=l 



m 



w 



2A1) 



~ 1 

m 



7,8 



Xi L 

VjlTkl — Vj2Tk2 



X] ■ 



\/2mw sin /3 



(1) 

7,8 \ ^2 



Xj ■ 



12 



a 



j2 f^W 



a/2 cos 13 m^. 
2 / 



(3) ( 



m 



Xi / 

mt 



k=l 



\plm-w sin (5 ^ 



Tkifr 



(3) f^t 



(38) 



A. Four exemplifying cases 



We are now ready to investigate the effects of the SUSY matchings on the b — > 577 
transition. We will show our results by plotting the ratios of SUSY versus SM decay rates 
for central values of the SM input parameters while varying the unknown SUSY parameters. 
We investigate four limiting cases which span the global features of the new amplitudes. 

First, we take M2 = = At = and fix tan/3 = 1 (case 1). In this approximation we 
have 



U 



V 



V2 



m 



1 -1 



XI, 2 



mw, 



(39) 



T 




~ 2 ~ 2 I 2 
m, ^ = m + m. 



The Wilson coefficients C7 g can be simply written as: 



(2x + ^)/«(x + ^)-^/S(x + z), 



where 



m 



X 



m 



2 ' 



(40) 



(41) 



(42) 

''W IT^w 

As remarked in ref. Cf^'^^(mvK) shows an exact cancellation in the supersymmetric 
limit, z —>■ and rriH — ^ mw- This is a consequence of the fact that any magnetic moment 
transition vanishes in exact supersymmetry |2^. Therefore non- vanishing contributions to 
the ^"^^ coefficients arise due to the presence of the soft breaking terms. 

We define 
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R 



'77 



and 



BR{B, ^ 77) 5M 



BR{B ^ X,7) SUSY 



(43) 



(44) 



where the SM decay rates are those given in eqs. (pOD -fpTf), obtained using the central values 



of the SM input parameters. 



Ryy J , 




200 

250 

m 300 




Fig. 4. Case 1. The allowed range for R^^ is shown as a function of rh and rriH (GeV), 
by constraining the BR{B —>■ Xs'~f)susY to vary within a ±30% from its SM value. 

In Fig. 3 we show R^^ as a function of the charged Higgs mass and the scalar soft 
breaking mass in a few hundred GeV range. In Fig. 4 the same range is spanned assuming 
the constraint 

0.7 <Rj< 1.3 . (45) 
14 



We see that R^^ as well is bound to vary in approximately the same range. The study of 
the ratio R^^/R^ in the same region shows deviations of at most ±4% from unity, which 
shows the strong correlation between the two decays. 

By releasing the constraint At = 0, while holding M = fi = and tan/3 = 1, we allow 
for a mass splitting of the stop eigenstates (case 2). This corresponds to having 



1 1' 



T 



The chargino contribution to Cy^g becomes: 



~ 2 ~ 2 I 2 



' Il,2 



+ m? ± Atfhmt . 



where Wi^2 = x + z ± Atmrrit/m 



w 



m 300 
250^ 

200^ 




(46) 



(47) 



Fig. 5. Case 2. R^-^ is shown as a function of m (GeV) and At for mn = 150 GeV. 

In Fig. 5 we plot i?^^ as a function of At and fh for fixed mn = 150 GeV, under the 
requirement that the lightest stop mass is always above 45 GeV. As we see, releasing the 
stop squark degeneracy while keeping charginos degenerate does not sizeably modify the 
features shown in case 1, and the same conclusions apply. 

Next we consider M = = = 0, and arbitrary tan (3 (case 3). The chargino component 
reduces to: 



X + z 
'4cos4/5 
z 



(1) 



4 cos^ (3 



(1) 



X + z 
2 cos2 (3 



2 cos2 f3 



^ A2) / x + z 



1 

2- 
1 

2' 



^ .(1) 



X + z 
2sin2/? 



(48) 
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In this case, the chargino degeneracy is hfted, while keeping the degeneracy in the squark 
sector. The chargino contribution becomes dependent on tan /3. On the other hand, as can 
be verified by means of eqs. (|30|)-(|35D, the tan/3 dependence of eq. ( ^Sf ) in the large tan/3 
limit is only logarithmic. As we will see later, for the SUSY amplitude to exhibit a stronger 
tan P dependence 7^ is required as well. 




Fig. 6. Case 3. The allowed range for R^^ is shown as a function of m and rriH (GeV) 



for tan/3 = 10, imposing the constraint of eq. (^51) on BR{B — > Xs'y)susY- 



300 m 
^250 




Fig. 7. Case 3. The allowed range for R^j is shown as a function of m (GeV) and tan/3 
for rriH = 150 GeV, imposing the constraint of eq. (ESf) on BR{B Xs'y)susY- 



In Figs. 6 and 7 we show as a function of different SUSY parameters the allowed range 
for R^^ once the constraint on BR{B — * Xs'y)susY in eq. P5| ) is imposed. R^^ is always 
bound to vary within ±30% from its SM expectation with high correlation to R^. 

Finally we consider the case for which M2, fi and At are different from zero and tan (3^1 
(case 4). In fact, the part of the chargino contribution which leads the large tan/3 behaviour 
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vanishes when either squarks or charginos are degenerate, as one can verify from the sim- 
phfied form of eq. (pS]) . 




Fig. 8. Case 4. The potential enhancement of R^^ is shown as a function of fh and nin 
(GeV) for tan/? = 15 (upper surface) and 10 (lower surface). 

An analytic approximation of the HulH^^ component of the chargino amplitude can be 
derived which shows explicitly its interesting features p3[. We assume the chargino mass 
matrix in eq. ( p6|) to be approximately diagonal: 

^ diag(M2, -/i) (49) 



This approximation holds effectively when |23| |M| — = 0[max(M|, /i^)] ^ and 
(M|, /i^) ^ m^. It is important to notice that these requirements, and therefore the 



approximation of eq. ( ^9]) is consistent with one of the eigenvalues, say being of the 
order of m^i/, while the other IM2I is much heavier. Thehe chargino mass matrix already 
being diagonal, the approximate mass eigenvalues are simply given by the absolute values of 
the parameters M2 and /i, and the two unitary rotations which "diagonalize" the chargino 
mass matrix can be written as: 

U ^ diag(sign[M2],-sign[/i]) , 

(50) 

V ^ 1 



In this approximation, the matrix T and the stop mass eigenstates are given by eq. (Pq). 
Using eq. ( ^Bj ) and eqs. (pUD-(^) we obtain a simple expression for the part of the 



chargino component relevant for large tan/?: 
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2 sin 2/? /i 



1 rut 




A3) 
J 7,8 




) 



(51) 



where rh^2 = l/^l is the hghtest chargino eigenvalue. Notice that the amphtude in eq. (|5TD 
depends on the signs of both fi and the trihnear soft breaking parameter At (changing the 
sign of the latter amounts to interchanging the two stop mass eigenvalues). One also verifies 
that the amplitude vanishes for either = or = as it should. 

As already mentioned, at variance with the case 3, the leading behaviour of the chargino 
amplitude is linear with tan f3. This may in general be the source of large deviations of the 
SUSY rates from the corresponding SM expectations. In Fig. 8 we show the ratio of the 
SUSY to SM branching ratio for — 77 as a function of mu and rh for tan (3 = 10 and 
15. In the example shown we have chosen /i = mw and At = 1.5. We see the potential large 
enhancements which arise for large tan/5 from this component of the chargino amplitude. 
On the other hand, imposing the constraint in eq. (|45|) allows only those regions of the lower 
surface for which R^^ varies approximately in the range 0.7 — 1.3. 

Globally, in the tested region of parameters the deviations of the Bg 77 decay rate 
from the SM expectations are confined to be within 10% from the corresponding deviations 
for the B Xg'y decay. 

At the next-to-leading order one may try to devise models that enhance the matchings 
of the O3-Q penguin operators (which are vanishing at the LO) keeping the Ot^s Wilson 
coefficients "under control". As unlikely as this may be, numerically it is anyhow difficult 
to expect drastic deviations from the SM predictions, due to the subleading role of the Os-q 
operators (analogous considerations apply to the electroweak penguin operators, which we 
have neglected in the LO analysis). 

We conclude that in order to disentangle new physics effects from a comparison of the 
two b ^ s radiative decays a precision below 10% is required both on the theoretical and 
experimental sides. Due to the smallness of the two-photon rates and to the theoretical 
uncertainties related to long-distance physics it shows as a truly challenging task. 
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TABLES 



as{mz) 0.118 

aem 1/129 

mz 91.19 GeV 

mw 80.33 GeV 

mt 175 GeV 

rub 4.8 GeV 

nic 1.4 GeV 

rus 0.150 GeV 

0.976 

\V*Vtb\ 4 X 10-2 

mfi^ 5.37 GeV 

fs. 0.2 GeV 

r{Bs) 4.09 X 10-1^ GeV 

DR{B — V,r;/, ) 10.1 x Kr- 



TABLE I. Values of the input parameters used in the numerical calculations. 
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